Hybrid superconducting quantum magnetometer 



F. GiazottcQ and F. Taddefl 

NEST Istituto Nanoscienze-CNR and Scuola Normale Superiore, 1-56127 Pisa, Italy 

A superconducting quantum magnetometer based on magnetic flux-driven modulation of the 
density of states of a proximized metallic nanowire is theoretically analyzed. With optimized ge- 
ometrical and material parameters transfer functions up to a few mV/$o and intrinsic flux noise 
~ 10 _9 $o/V / Hz below 1 K are achievable. The opportunity to access single-spin detection joined 
with limited dissipation (of the order of ~ 10 -14 W) make this magnetometer interesting for the 
investigation of the switching dynamics of molecules or individual magnetic nanoparticles. 



I. INTRODUCTION 

The superconducting quantum interference device^ 
(SQUID) is recognized as the most sensitive magnetic- 
flux detector ever realized, and combines the physical 
phenomena of Josephson effect 2 and flux quantization^ 
to operate. SQUIDs are nowadays exploited in a vari- 
ety of physical measurements^ with applications span- 
ning, for instance, from pure science to medicine and 
biolog}BEl Recently, the interest in the development of 
nanoscale SQUIDs 7 9 has been motivated by the oppor- 
tunity to exploit these sensors for the investigation of 
the magnetic properties of isolated dipole a 10 " 14 ', with the 
ultimate goal to detect one single atomic spin, i.e., one 
Bohr magneton. 

Here we theoretically analyze a hybrid superconduct- 
ing interferometer which exploits the phase dependence 
of the density of states (DOS) of a proximized metal- 
lic nanowire to achieve high sensitivity to magnetic flux. 
The operation of a prototype structure based on this 
principle, the SQUIPT 15 , has been recently reported. 
We show that with a careful design transfer functions 
as large as a few mV/$o and intrinsic flux noise ~ 
10 _9< I ) o/a/Hz can be achieved below 1 K. Limited dissi- 
pation joined with the opportunity to access single-spin 
detection make this structure attractive for the investi- 
gation of the switching dynamics of individual magnetic 
nanoparticles. 

The paper is organized as follows. The model of the hy- 
brid superconducting magnetometer is presented in Sec. 
II. The Josephson and quasiparticle current are calcu- 
lated in Sees. Ill and IV, respectively. The flux resolu- 
tion and device performance are finally presented in Sec. 
V, where we address briefly the feasibility of this struc- 
ture as a single-spin detector. Sec. VI is devoted to the 
Conclusions. 



II. MODEL 

The interferometer [sketched in Fig. [TJa)] consists of 
a diffusive normal metal (N) wire of length L in good 
electric contact (i.e., ideal interface transmissivity) with 
two superconducting electrodes (Si) which define a ring. 
We assume the wire transverse dimensions to be much 
smaller than L so that it can be considered as quasi-one- 
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FIG. 1. (Color online) (a) Scheme of the device. L is the 
wire length whereas w is the width of the superconducting 
tunnel junction (S2) coupled to the middle of the N region. 
ip is the macroscopic quantum phase difference in Si, while 
$ is the magnetic flux threading the loop. Furthermore, Rt 
is the tunnel junction normal-state resistance and hias is the 
current flowing through the structure. 



dimensional. The contact with Si induces superconduct- 
ing correlations in N through proximity effect! 16 " 20 ! w hich 
is responsible for the modification of the wire DO!9^4 For 
lower-transparency NSi interfaces the proximity effect in 
the wire will be reduced thus weakening the effects de- 
scribed below. In addition, a superconducting junction 
(S2) of width w and normal-state resistance R t is tunnel- 
coupled to the middle of the N region. The loop geome- 
try allows to change the phase difference (cp) across the 
NSi boundaries through the application of a magnetic 
field which modifies the wi re DO and the transport 
through the tunnel junctiori 15 l 23 i 

The proximity effect in the wire can be described with 
the quasiclassical Usadel equations^. The short-junction 
limit (i.e., for Ai <C hD/L 2 = E T h, where Ai is the or- 
der parameter in Si, D is the wire diffusion constant, and 
Eth is the Thouless energy) will be considered in the fol- 
lowing, since in such regime the Usadel equations allow 
an analytic expression for the wire DOS 24 thus simplify- 
ing the device analysis. In addition, the interferometer 
performance is optimized in this limit as proximity effect 
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in the wire is maximize d 22 * 24 *. Assuming a step-function 
form for the order parameter Ar^Sl, i.e., constant in Si 
and zero in the N wire, the wire DOS normalized to the 
DOS at the Fermi level in the absence of proximity effect 
is given by 2 ^ 



N N (x, s, T, ip) = Re {cosh [0(x, e, T, <p)]} , (1) 



where 



= arcosh[a(e, tp, T)cosh[2xarcosh[/3(£, <p, T)]]) , (2) 



and 



V[s 2 -A2(T)cos2(<p/2)] (3) 



[e 2 - Af(T)cos 2 (y>/2)]/[e 2 - Af(T)]. (4) 



In the above expressions, e is the energy relative to the 
chemical potential of the superconductors, T is the tem- 
perature, and x G [—L/2, L/2] is the spatial coordinate 
along the wire. Nn exhibits a minigap (e g ) 



7 (^) = A 1 (T)|cos(^/2)| 



for |e| < e g whose amplitude depends on cp and is con- 
stant along the wire. In particular, e g = Ai for cp = and 
decreases by increasing vanishing at = 7r. Finally, 
by neglecting the ring inductance the phase difference 
becomes ip = 27r$/<I>o, where is the total flux through 
the loop area, and $o = 2.067 x 1(T 15 Wb is the flux 
quantum. 



III. JOSEPHSON CURRENT 

At equilibrium a current through the system I eq can 
flow thanks to a direct Josephson coupling between 
the superconducting electrode S2 [with order parameter 
equal to A2(T)e* Xs where xs is the macroscopic phase] 
and the proximized N wire. We assume the BCS tem- 
perature dependence for Ai )2 (L) with critical temper- 
ature T c i 5 2 = A? 2 /(1-764&b) where A? 2 is the zero- 
temperature order parameter in S 1,2, and is the Boltz- 
mann constant. Since the junction NS2 is extended in the 
^-direction one can calculate, in the tu nneling limit, I eq 
using the quasiclassical approac h 1 16 * 27128 1 as the following 
integral 
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I rw/2 
— — / dx 

XjR t J-w/2 



(6) 
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where cr 3 is the third Pauli matrix, [•, •] represent the 
commutator, and e is electron charge. Furthermore, 
Gj[(e,x,(p) is the position-dependent retarded Green's 




FIG. 2. (Color online) Josephson current vs flux (I eq — <£) 
characteristics calculated for a few values of temperature T 
assuming A J = O.lErh and A J = 4A°. Thick lines represent 
the supercurrent calculated using the Ambegaokar-Baratoff 
formula Eq. (11 ), whereas the thin lines are the supercurrent 



calculated for an extended tunnel junction with width w 
L/2. 



(5) function on the N wire and G^ 2 (e) is the retarded Green's 



function of the electrode S2. They are defined as follows: 



and 



G%(e,x,tp) 



cosh 9 sinh 0e lx 
sinh#e _ * x — cosh# 



e A 2 e ixs 



(7) 



(8) 



where is defined in Eq. ([2| and 

X = — arctan[7(e, T)tanh[2xarcosh[/3(£, (p, T)]]] , (9) 
with 

7 = sj[e* - A?(T)cos2( ¥ ,/2)]/[A 1 (T)cos( V /2)]. (10) 

Note that the supercurrent I eq depends on the phase xs 
of the order parameter in the electrode S2. We are inter- 
ested in the critical current that we determine by fixing 
Xs such that it gives the maximum supercurrent for a 
given value of cp. In Fig.|2]the equilibrium critical current 
is plotted, for different values of temperature, as a func- 
tion of cp in units of A^/ieRt), assuming Aj = 0.1^, 
A^ = 4A^ and w = L/2. For the sake of compari- 
son we also plot the supercurrent calculated through the 
Ambegaokar-Baratoff 29 formula, relative to a point-like 
NS2 junction: 



rAB 



Tre g ( l p)A 2 (T)k B T 
eR t 



(11) 



Z=0,±1,- 
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where uoi = nk^T(2l + 1). Interestingly, for our choice of 
parameters the difference between I eq and is hardly 
noticeable: the lateral spatial extension of the NS2 junc- 
tion along x plays a marginal role. As a matter of fact 
the supercurrent turns out to be negligible in the exper- 
iments reported in Refs. [T5landl37l The reason for such 
Josephson current suppression is presently unclear. One 
possibility is that the system is brought out of equilib- 
rium by voltage fluctuations which might originates from 
the measuring circuit. Such voltage fluctuations will drop 
mostly across the tunneling barrier, being the most resis- 
tive component of the system. As a result, the Joseph- 
son current will oscillates at high frequency, hindering 
the possibility of detection. This fact is fortunate, since 
a supercurrent might prevent a correct voltage read-out 
of the device, which is not observed experimentally. 



IV. QUASIPARTICLE CURRENT 

The current through the tunnel junction biased at volt- 
age V is therefore dominated by the quasiparticles, and 
can be written as^ 



I rw/2 r 

1= — 5- / dx / deN N (x,e,<p)N S 2(e)F(e,e), 
ewR t J_ w/2 J 



where 



N S2 (s,T) 



\A 2 -A 2 (27 



(12) 



:6[£ 2 -A 2 (T) 2 ] (13) 



is the normalized DOS of the S2 electrode, e — e — eV, 
Q(y) is the Heaviside step function, F(e,e) = [fo(e) — 
fo(e)], and f (e) = [1 + exp^/Zc^T)] -1 is the Fermi- 
Dirac energy distribution. In the following we set A2 = 
200 /ieV and A? = 4A§ = 800 /ieV as representative val- 
ues for a structure exploitin g alum inum (Al) and vana- 
dium (V) as superconductor J 3Q I 31 I, respectively, w = L/2 
and R t = 5 Mfi. 

Figure 3 shows the interferometer current vs voltage 
(I-V) characteristics calculated at T = 0.1T C 2 for dif- 
ferent values of the applied flux <E 32 . In particular, for 
$ = the I — V characteristic resembles that typi- 
cal of a superconductor-insulator-superconductor junc- 
tion (i.e., S1IS2 where I denotes an insulator) where the 
minigap in the wire is maximized [i.e., e g = Ai(T)l, 
and the onset for large quasiparticle current occurs atP 
V = [Ai(T) + A 2 (T)]/e. For $ = $ /2 the charac- 
teristic is similar to that of a normal metal-insulator- 
superconductor junction (i.e., NIS2) with e g suppressed. 
The curves show a peak at V = |Ai(T) - A 2 (T)|/e which 
corresponds to the singularity appearing in the tunnel- 
ing characteristic between different superconductors 3 -^. In 
a current-biased setup the interferometer operates as a 
flux-to-voltage transducer providing a voltage response 
V(<f>) that depends on the bias current ifo as fed through 
the tunnel junction [see Fig. 3]. For any I^ as , V(3>) is 
determined by solving the equation Iuas — 1 = 0. 
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FIG. 3. (Color online) Interferometer quasiparticle current 
vs voltage (7 — V) characteristics calculated for a few values 
of <3> at T = 0.1T C 2. T C 2 is the critical temperature of S2, 
Iuas is the current flowing through the device and V($) is 
the resulting voltage modulation. 
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FIG. 4. (Color online) (a) V vs $ calculated for several bias 
currents hias at T = 0.1T C 2. (b) T vs $ calculated for the 
same hias values and T as in (a). 



Figure 4(a) shows at T = 0.1T c2 calculated for 

several Ibias values. V(&) turns out to be maximized at 
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FIG. 5. (Color online) (a) V vs <3> calculated for a few tem- 
peratures at hias = 1.0 A° /eRt. (b) T vs <E> calculated for the 
same T values and hias as in (a). 



the lower bias currents where the voltage swing obtains 
values as large as ^A^/e, whereas it is gradually reduced 
by increasing hias- The interferometer performance is 
thus improved at low hias- 

An important figure of merit of the interferometer is 
represented by the flux-to-voltage transfer function 1 



dV 
<9$ 



(14) 



which is shown in Fig. 4(b) for the same hias values as in 
panel (a). In particular, T as large as ~ 12.5A2(e<E>o) _1 
can be obtained at the lowest currents, whereas it is grad- 
ually suppressed at higher hias- 

The role of the temperature is shown in Fig. 5(a) 
which displays V(<&) calculated for several T values at 
/ = 1.0A2/(ei?t). An increase in T leads to a reduction 
of V($) as well as to a suppression and smearing of the 
voltage swing. This directly reflects on the transfer func- 
tion, as displayed in Fig. 5(b). We note that even at 
T = T C5 2, i.e., when S2 is driven into the normal state, 
T as large as c± S.GA^e^o) -1 can be achieved. It fol- 
lows that voltage swings up to 0.8 mV and T as large as 



2.5 mV/$o can be achieved with the suggested materials 
combination for T < 1 K. 



V. NOISE AND DEVICE PERFORMANCE 

We now turn on discussing the noise properties of the 
interferometer. In the actual current-biased setup an im- 
portant quantity is represented by the voltage noise spec- 
tral density (Sy) denned as 



R 2 d Si, 



(15) 



where Rd = dV/dl is the tunnel junction dynamic resis- 
tance, and Si is the current noise spectral density (shot 
noise) given 



2 rw/2 p 

Si = —BT dx dsN N (x,e,$)N S2 (e)M(s,s), 

wRt J- w /2 J 



(16) 



where 



M{e,e) = / (e)[l - /„(e)] + /„(e)[l - f (s)]. (17) 

The intrinsic flux noise per unit bandwidth of the inter- 
ferometer ($ nj ) is related to the voltage noise spectral 
density aa^ 



1^)1' 



(18) 



Note that <£ ns oc \fRt, as Sy oc Rt and F($>) is indepen- 
dent of tunnel junction resistance. 

Figure [6^a) shows <£ ns versus hias for several flux val- 
ues at T = 300 mK. <£ ns is a non-monotonic function of 
hias with a minimum which depends, for each <3>, on the 
bias current. In particular, an increase in <$> leads to a 
general reduction of $ ns at low hias, while its minimum 
moves toward lower bias current. We stress that <£ ns as 
low as 10 -9 ^o/a/Hz or better can be achieved at this 
temperature in the ~ 10... 80 pA range for suitable val- 
ues of <£. This good flux sensitivity stems from the low 
shot noise Si (which is peculiar to all-superconducting 
tunnel junctions) together with a small Rd at the biasing 
point and large J 7 ^). 

The temperature dependence is displayed in Fig. |6|b) 
where <£ ns vs hias is plotted for different T values at 

= 0.3^0- Notably, the minimum of <£ ns turns out 
to be quite insensitive to the temperature up to ~ 900 
mK. Then, higher T yields to a reduction of the current 
window suitable for high flux sensitivity and to an over- 
all enhancement of <£ ns . Furthermore, for T > T C 2 [see 
the line corresponding to T = 1.5 K in Fig. 6(b)] & ns 
is significantly degraded in the whole hias range. This 
emphasizes the effectiveness of a superconducting tunnel 
probe for a drastic suppression of <£ ns . 

The impact of dissipation P = VI is displayed in 
Fig. [6jc) which shows P vs hias f° r different values at 
T = 0.3 K. P can largely change by varying and hias 
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FIG. 6. (Color online) (a) <3> n s vs hias calculated for differ- 
ent values at T = 0.3 K. (b) <£ ns vs hias for $ = O.3<3> 
calculated at different temperatures, (c) P vs hias calculated 
for different <E> values at T — 0.3 K. In all calculations we set 
A° = 200 fj,eV, A? = 800 /xeV, and ft = 5 MQ. 



Ag, A 1 /E Th ~ 0.3, and 7 C ° ~ 318 /iA. By choosing, for 
instance, a circular washer geometr^ with 2r = 150 nm 
as internal diameter and external radius ^ > r we get 
£g = ~ 0-19 pH, where /iq is the vacuum perme- 

ability, so that 27r/ c £ G /$o ~ 0.18. Condition ii) can be 
fulfilled as well by choosing a suitable washer thickness 
and 7Z. 

As Cq has to be kept small to satisfy condition i) it fol- 
lows that the present structure could be suitable for the 
measurement of the magnetic properties of small isolated 
samples. In this context, the magnetometer sensitivity 
(<S n ) to an isolated magnetic dipole place d at the center 
of the loop is approximtely given b y 1Q l 12 l 



(19) 



where \±b is the Bohr magneton. With our choice for 
r and by coupling the device to a cryogenic voltage 
preamplifier^ (which we assume dominates the volt- 

^ 0.1 uV/a/Hz yields a to- 

m^Z)\ ~ 40 ^o/VRz, lead- 
ing to S n « 1 atomic spin/V^Hz below 1 K. Further- 
more, the best achievable energy resolution 1 would be 



age noise) with y^S 
tal flux noise <&t°* 
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as well. In particular, in the ~ 10... 80 pA current range, 
P can vary from a few fW to some tens of fW. Such a 
small power has the additional advantage to prevent sub- 
stantial electron heating in the N wire^H By contrast, in 
conventional SQUIDs dissipation is typically from two to 
five orders of magnitude large ]™. As P oc Rt 1 , dissi- 
pation can be tailored by choosing a proper value of the 
tunnel junction resistance. 

For a correct operation of the interferometer the two 
following conditions should be fulfilled: i) 2ttI^Cg/^o i$ 
iS (where I® is the zero-temperature critical current 
of the SiNSi Josephson junction, and Cq is th e loop 
geometric inductance), and ii) C^ 1 <C £^ 15 l 23 l [where 
C^ ljN ~ HRs 1 ,n/^^i is the kinetic inductance^, and 
Rs!,n is the normal-state resistance of Si(N)]. Condi- 
tion i), where 7° = OMiTA^/eR^ ensures to avoid 
magnetic hysteresis whereas ii), which is equivalent to 
Rs 1 <C Rn i ensures that the phase difference set by 

drops entirely at the wire ends thus allowing a full 
modulation of its DOS. As an additional set of pa- 
rameters we choose a silver (Ag) wire with L = 80 
nm, cross section A = 30 x 10 nm 2 , and D = 0.02 
m 2 s _1 which yield Rn = Lj{Av F e 2 D) ~ 5.2^, where 
i/ F = lx 10 47 J^m -3 is the DOS at the Fermi level in 



VI. CONCLUSIONS 

In summary, we have theoretically investigated a hy- 
brid superconducting magnetometer whose operation is 
based on magnetic flux-driven modulation of the density 
of states of a proximized metallic nanowire. In partic- 
ular, we have shown that with suitable geometrical and 
material parameters the interferometer can provide large 
transfer functions (i.e., of the order of a few mV/$o) and 
intrinsic flux noise down to a few n^o/V^Hz below IK. 
Furthermore, joined with limited power dissipation, the 
structure has the potential for the realization of sensi- 
tive magnetometers for the investigation of the switching 
dynamics of small spin populations. 
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